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This article extends the work of Bernabeu and Espinoza by examining the CP-violation reach 
of a ^^°Dy electron capture beam through the variation of the two Lorentz boosts, the number of 
useful electron capture decays, the relative run time of each boost and the number of atmospheric 
backgrounds. The neutrinos are assumed to be sourced at CERN with an upgraded SPS and are 
directed towards a 440 kton Water Cerenkov detector located at the Canfranc laboratory. Two 
' large 'CP-coverage' choices for the boost pairings are found; a (5-symmetrical coverage for (71,72) 

= (280, 160) and an 5-asymmetric coverage for (71,72) = (440,150). With a nominal useful decay 
rate of Nious = 10^* ions per year, the (5-symmetric setup can rule out CP-conservation down to 
sin^ 2^13 = 3 • 10~*. To reach sin^ 2S13 = 1 ■ 10^'^ for both 5 < and (5 > requires a useful decay 
rate of Mons ~ 6 • 10^^ ions per year. 
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I. INTRODUCTION 

The results of a series of ground breaking neutrino experiments indicate that neutrinos are both massive and that 
pi • they mix amongst themselves. Data from atmospheric [1, 2], solar [3-5], reactor [6-8] and long-baseline accelerator [9, 
O \ 10] neutrino experiments can be accommodated by two approximate 2-neutrino mixing schemes, each parameterised 
i-S^ . by a mass-squared splitting and a mixing angle. 

The goal of the future experimental program in neutrino oscillations is to determine how these 'solar' and 'atmo- 
. spheric' sectors combine into an overall 3-neutrino mixing scheme. This requires a search for sub-dominant appearance 
J> ' events; a signal whose strength is controlled by a third angle, 6*13, that is currently only bounded from above. Up 
, to three physical phases can be incorporated into a 3-neutrino mixing scheme; however, only the effect of the Dirac 
■ phase, 5, can be seen in an oscillation experiment as CP-violation. If the third mixing is zero then the two 2-neutrino 
\ mixing schemes are disjoint and there is no discrepancy between the neutrino and anti-neutrino mixing matrices 
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Defining Am^j = — m| , the current best fit values for the oscillation parameters are [11] 



O ■ Am 



jAm^il = 2.4 X 10"^ eV^ and sin^ 6I23 = 0.50 ; 
Am^i = 7.65 X 10"^ eV^ and sin^ 812 = 0.304 . 



The third mixing angle, ^13, is constrained to be [11] 



sin^ 6*13 < 0.040 (0.056) at 2a (3cr) ; (1.1) 

although some collaborations report a hint at low significance from a combined analysis of atmospheric, solar and 
long-baseline reactor neutrino data [12, 13]. 

Running and near future neutrino oscillation experiments [10, 14-16] will be the first to probe ^13 below the current 
limit and possibly confirm the reported hint for 6*13 = 0. If ^13 is beyond the reach of near future experiments, 
intense sources of neutrinos will be necessary: next generation SuperBeams [17], Neutrino Factories [18, 19] and Beta 
Beams [20, 21]. The aspiration of future long baseline neutrino oscillation experiments is to determine the third 
mixing angle, 6*13, whether there is CP- violation in the lepton sector; and to determine the sign of Am|]^. This is not 
straightforward; the analysis of the appearance channels is capable of returning up to 8 equally allowable fits to the 
data [22-26]. The challenge is therefore successful resolution of these degeneracies and a simultaneous push for the 
best sensitivity. 

A Beta Beam is a source of intense and coUimated beams of electron neutrinos or anti-neutrinos originating from the 
decay of boosted radioactive ions. The proposal is an integrated facility combining individual production, acceleration 
and storage components. Ion production could make use of the Isotope Separation Online (ISOL) technique [27] within 
the EURISOL at CERN; production rings [28] to minimise losses of parent nuclides in a beam dump; or a two target 
system using an intense beam of deuterons [29] . The acceleration would make use of existing or potential upgrades to 
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current facilities such as the CERN accelerator chain, the Fermilab Main Injector and Tevatron; and maybe DESY. A 
storage ring whose straight sections point towards the far detectors will be the main infrastructure addition required. 

It is advantageous to produce a clean, coUimated beam by this technique since a single boost is all that is necessary 
for a laboratory neutrino flux to span an energy range up to several GeV. However, one is not free to choose the 
laboratory neutrino spectrum; the boost determines the maximum neutrino in the laboratory frame which in turn 
determines the un-oscillated event rate at all energies. One way to alleviate this problem is to use a source of mono- 
energetic neutrinos so that the boost of the ion can be altered to freely choose the laboratory frame energy of the 
neutrinos. Clearly, such an approach needs at least two boosts as both ^13 and S are unknown. 

The use of electron capture decays to source mono-energetic neutrino fluxes has taken on two guises. The first 
paper [30] (and a subsequent study [31]) proposed the use of ions with Q- values less than twice the electron mass, 
hence positron decay is kinematically forbidden. These are not realistic, however, since the ion half-lives are too long 
and the boosts required can only be achieved with the LHC The idea of using electron capture decays was in fact 
first suggested at the EURISOL meeting in January 2005 [32] with focus on higher Q-valued ions. This was later 
expanded into studies for the CERN-Prejus and CERN-Canfranc baselines [33, 34]. For ions with Q- value greater than 
twice the electron mass, the electron capture channel competes with a positron decay 'background'; the branching 
ratio for electron capture channels drops sharply as the Q-value increases. With Q-values ~ 3-4 MeV and a 1 TeV 
accelerator, one would be able to placx; a mono-energetic neutrino beam on first oscillation maximum for baselines 
up to 1000 km, making the idea an attractive one in Europe with the source based at CERN. However, the positron 
decay background would make the idea unattractive since the number of useful electron capture decays would not be 
high enough for a competitive physics reach. It could be argued that the positron decay neutrinos could be used to 
boost sensitivity by providing coverage of the second oscillation maxima with concentration of the electron capture 
neutrinos on or around first oscillation maximum. This is the hybrid Beta Beam strategy explored in detail in [35]. A 
number of nuclides that decay quickly through large Gamow- Teller resonances have been discovered in the Gd region 
of the Segrc chart, see for example [37], that opens the possibility of a higher Q-value electron capture beam. These 
ions source dominant electron capture channels at higher Q-values. One of these ions (^■''°Dy) forms the basis of this 
paper as it did [33, 34] where it was used to source the CERN-Frejus and CERN-Canfranc baselines with a large 
Water Cerenkov as the far detector. 

In this article, the work of Bernabeu and Espinoza [34] is extended to include the effect of the choice of Lorentz 
boost pairing, the relative run times at each boost, the number of useful ion decays, and the number of atmospheric 
neutrino events. In Sec. II, the electron capture beam concept is reviewed with some discussion on the choice of ions 
and some of the technological difficulties not present in a standard Beta Beam. In Sec. Ill a variant on the intrinsic 
degeneracy, which I refer to as the energy degeneracy, is introduced which results from the combination of events 
at two energies for a neutrino or anti-neutrino run only. The analyses carried out are summarised in Sec. IV; the 
variation of the boosts is performed in Sec. V; and the number of useful decays and atmospheric backgrounds are 
included in Sec. VI. To finish, the use of mono-energetic neutrinos from bound beta decay is discussed in Sec. VII, 
and the conclusions are drawn in Sec. VIII. 



II. THE ELECTRON CAPTURE BEAM CONCEPT 



Electron capture is a decay channel available to proton-rich nuclei and it competes with positron decay depending 
on the energy available. For Q < 2me, positron decay is kinematically forbidden and so electron capture decays form 
the entire phase space. For electron capture decays, the rate is proportional to the from the two body decay and 
from the square of the electron orbital wave-function; where Q is the Q-value of the decay and Z is the proton 
number. The relative rate with respect to positron decay for Q > 2me is approximately given by 

where the matrix elements for the channels are assumed to be identical, a is the fine structure constant, and the 
leading is taken for Tp. For electron capture with Q-value Qec and boost 7, the neutrino flux in the laboratory 
frame at baseline L from source is given by [30, 33, 34] 

'^^ ^'""^ 7' 5(27Qec - E,) = (5(27Qec - E,) . (2.2) 



dQdE,, 7rL2 

Trivially, the baselines available to electron capture machines are dependent on the Q-value of the ion and the 
maximum boost allowed from the acceleration. 

A collection of ions in the Gd region of the Segre chart, suitable for European baselines, are presented in Tab. I. 
i50j3y T^g^g chosen in [33, 34] as the resonance did not have a width and because of its relatively low Q-value. From 
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TABLE I: Candidate electron capture beam ions. Based on a similar table from [37]. 



Eq. 2.2, the neutrino flux is proportional to the square of the boost for a fixed baschne and number of useful decays. 
Lower Q- values require higher 7 to achieve the same laboratory neutrino energies, and so lower Q- values are favoured 
in this sense, provided the accelerator is capable of achieving the boosts. The discovery of ions far from the stability 
line that decay through a giant Gamow- Teller resonance opened up the possibility of using 1 TcV machines such 
as an upgraded Super Proton Synchrotron (SPS) [36] (envisaged in some LHC-upgrade scenarios), to source the 
intermediate baselines in Europe. As mentioned previously, [33, 34] considered this probability for the CERN-Frejus 
(130 km) and CERN-Canfranc (650 km) baselines with ^^''Dy as the chosen ion. Since only a neutrino flux is available 
in an electron capture setup, the approach was to exploit the different energy dependence of the CP-odd and CP-even 
properties of the appearance probability (Eq. 3.2) by running at two different boosts. The purpose of those studies 
was not to identify Lorentz boosts of particular interest, rather to demonstrate the phenomenological feasibility of 
the idea: the use of multiple boosts to extract the energy dependence of the event spectra. The goal of this article is 
to examine the CP-violation reach as a function of the boosts, number of useful decays and atmospheric event rate. 

For electron capture and Beta Beam setups, the parent ions will be accelerated in the existing or upgraded CERN 
infrastructure before accumulated and stored in a ring whose long, straight sections source the neutrino flux. For 
electron capture machines with ion boost 7, the mono-energetic neutrino flux at a detector distance L from the source 
is given by Eq. 2.2. Since the neutrino energy in the laboratory frame is given by E^, = 2jQeC: for fixed baseline and 
number of useful decays, the ions with lower Q-values result in larger fluxes (because of the larger boost needed to 
achieve a given energy). The choice of ion is thus a balance between Q-value and the available acceleration. From the 
ions presented in Tab. I, ^*^Dy and ^^*Er have very similar characteristics; the bulk of all decays being an electron 
capture, half-lives ~ 3 minutes and neutrino energies of ~ 2 MeV in the rest frame. The lower Q-value ion, ^^"^Dy, will 
be taken in this paper for consistency with [34] and the higher 7's required compensate for the lower electron capture 
branching ratio. The remaining 36 % for ^^''Dy is a-decay and so does not source a primary neutrino background; the 
daughter has 100 % Q-dccay with a 74 year half-life. The other ions in Tab. I have undcsircable intensity and Q-value 
combinations. These higher Q-value ions also have a large positron decay background and are more suitable for the 
hybrid machine introduced in [35]. For a ^^"^Dy ion, the maximum boost attainable with an upgraded 1 TeV SPS [36] 
is 7max = 440^ corresponding to a laboratory frame energy of = 1.23 GeV. For a electron capture machine with 
no backgrounds the choice of detector does not depend on the energy reconstruction capabilities as the event energy 
is determined by the boost. In reality, some reconstruction might be necessary for the appearence events to reduce 
the atmospheric background, especially if low production rates or issues with the acceleration force a large duty cycle. 
This point will be discussed further in Sec. VI. Since the maximum laboratory frame energy is 1.23 GeV, baselines 
in excess of CERN-Canfranc are unrealistic propositions. Therefore matter effects are small, but not negligible, and 
are not enough to achieve competitive sensitivities to the mass hierarchy. Electron capture machines of this type are 
therefore 'CP- violation machines'. 

Although an electron capture machine is a variant on a Beta Beam, there are number of further technological 
challenges that need to be overcome for such a facility to be realised; 

1. Candidate ions have large proton numbers which bring problems with space charge in the early part of the 
accelerator chain. 

2. Electrons need to be left on the ions for the electron capture channel to be available. 

3. Decays rates are slow; typically several minutes. The optimal half-life for a Beta Beam type machine is 0{1 
sec). 



The maximum boost for a fully stripped ^^'^Dy ion is in fact 7 = 468; however, it is necessary to leave several electrons bound to the 
nucleus to source the electron capture decay. If one leaves the 2 K-shell electrons and 2 in the L-shell, this reduces the maximum boost 

to 440 
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All these issues conspire to make the desired useful rate of 10"'^^ ion decays per year extremely hard. To date, no 
extensive R&D has been carried out on the possible yearly rates for the rare-earth nuclei sought for electron capture 
machines. However, code adapted from the baseline Beta Beam study has been used to simulate estimate rates with 
the same vacuum and accelerator conditions [38]. Assuming a round Gaussian beam of ions, the self-field incoherent 
tune shift was found to be large. Increased intensities in the acceleration only worsen the space charge problem. This 
problem can be tempered by leaving extra electrons on the ions, thus reducing the charge; however, such a strategy 
reduces the maximum boost attainable. This is not necessarily bad, as the results of the simulation in Sec. V indicate 
that the maximum boost 7 = 440 may not be necessary. A possible strategy could be to leave electrons bound to the 
ion such that the maximum boost is the sought boost. For instance, if one requires a boost of 7 = 150, 44 electrons 
can be left bound to the Dy nucleus. The tune-shift has a charge-squared dependence, so we decrease its effect by a 
factor of 9 with this action. 

This, however, brings with it a new problem. The vacuum conditions of the accelerator chain and decay ring are 
such that the probability of electrons to be stripped from the ion is significant. Since this would alter the charge-to- 
mass ratio of the ion, the magnetic configuration of the accelerator or decay ring no longer matches that required and 
the ion is lost. These losses can be thought of as an extra decay channel. The annual rate of neutrinos needs to be 
modified [38]: 

^ ^ IjJ Aec/7 _ gmT,ep(AEc/7-KAvac)'\ rp^^^ 

^rep Aec/7 + A vac ^ ^ 

Here, Aec and Avac are the decay constants for the electron capture decay and vacuum losses respectively; Tj-op is the 
repetition period for fills in the decay ring; Jin the total number of ions injected into the decay ring for each fill; I 
the livetime (the fraction of the decay ring length that sources the neutrino beam); m is the number of merges in the 
decay ring and T^un the length of the experimental run in seconds. It is reasonable to assume 1/Avac ~ 60 sees [38]. 

Ultimately though, the main problem in sourcing a competitive number of useful decays is the half-lives of the 
ions. Sourcing a beam of mono-energetic neutrinos in this manner requires the use of ions with a fast Gamow- Teller 
resonance. This requirement is very restrictive in terms of the choice of ions. Out of the candidates identified, the 
fastest decays have weak EC-transitions; those with dominant resonances have much longer half-lives 3 mins). To 
realise an electron capture machine therefore will require a thorough search of the Gd region of the Segre chart for 
an ion possessing a suitably fast resonant electron capture decay. 



(2.3) 



III. THE 'ENERGY' DEGENERACY 



It is well known that the analysis of data from a future long baseline facility suffers from the problem of degen- 
eracies [22- 26]; the asymmetry between neutrino and anti-neutrino probabilities, the unknown sign of Arrt^i, and the 
unknown octant of O23 can all lead to multiple fits to experimental data. For binned data, the number of neutrino 
(anti-neutrino) events in the ith neutrino (anti-neutrino) energy bin for the pair (^13 , S) is given by 



Ni{hz, ~5)=NTt / e(^,)a,^(p^)(^,) P± (i;,, ^13, ~5) $.,(p,)(S,) dE^ , (3.1) 

JEi 

where A/r is the number of targets in the detector, t is the time of data taking, e{E^) is the detector efficiency, 
(t{E^) is the interaction cross section, ^{E^) is the beam spectrum and AE is the bin width. Using the shorthand 
Aji = Am.'jJ{2E), the oscillation probability Pi^^^i,^ = Pf,^ can be expanded in the small parameters ^13, A21/A31, 
A21/A and A21L [39], 

P± (^-13, ~5) = sin^ 2^-13 sin^ 0,, (^^) ' sin^ 

^^A2iA3i . }al\ . (B:^L\ ( - A31L 
+ J — T- -TT- sm Z, cos ±d — 



A \ 2 J V 2 

-f- cos^ ^23 sin^ 2^12 (^^) ^^'''(t^) ' (^-^^ 

where J = cos ^13 sin 2^12 sin 2^23 sin 2^13, the ± corresponds to neutrinos/anti- neutrinos and Sip = A31. 

Here we are using A = \/2Gpnc^{L) (the constant density approximation for the index of refraction) where 
fie = ^1 L ne{L')dL' is the average electron density and ne{L) is the electron density along the baseline. 



5 




FIG. 1: Equi-probability curves for the CERN-Canfranc baseUne (650 km) for sin6lf3 = 0.05 and 5^' = 50°. The other 
oscillation parameters have been set to their current central values [11]. The black curves use = 1.3 GeV and the red curves 
use Ev = 1.8 GeV. Solid lines are for neutrinos, dashed lines are for anti-neutrinos. 



Labelling A''"'" to be the number of neutrino events, and A'^ to be the number of anti-neutrino events; 4 pairs of 
equations can be solved 

iV±(^i3,J,|Am^i|,023) =iV±(0i3,(5,|Am2^|,023) , (3.3) 
A^±(^i3,^,|Am2i|,023) -iV±(0i3,(5,-|Am2,|,023) , (3.4) 
N^{0uJ,\Aml,l923) = N^{9i3,S,\Aml^l90" ~ 023) , (3.5) 
N^ie,3j,\Aml,ie23) ^ N^{0u,S,-\Aml,l90'' - 023) ; (3.6) 

which in general lead to 8 solutions that can fit the data. In order, the first set of equations results in the 'intrinsic 
clone'; the second returns the 'hierarchy clones'; the third gives the 'octant' clones; and the fourth equation allows for 
'mixed' clones. Solutions to the first set of equations are depicted graphically in Fig. 1 for sin^ 29i3 = 10~^, i5 = 50°, 
023 = 45° and L = 650 km. The black lines correspond to E = 1.3 GeV (first oscillation maximum) and the red lines 
are for E = 1.8 GeV. Solid lines are for neutrinos and dashed lines are for anti- neutrinos. 

From Fig. 1, there is another type of degeneracy not discussed in the literature. Consider the first set of equations 
above but now also for energies Ei and E2 : 



N^J9i3,~S,\^ml,\, 623) 



N^J9,3,5,\Aml,\,923) 



(3.7) 



Attempting to find a solution to these equations will return only the true solution since the location of the intrinsic 
clone is energy dependent. This is the main strategy in resolving the intrinsic degeneracy. Suppose, one does not 
have an anti- neutrino channel. This set reduces from 4 to 2 equations, and in general will possess a clone solution. 
This can been seen as the intersection of the red and black solid curves at {9i3,6) — (4.8°, —140°). The location 
of this degeneracy will be different for neutrino and anti- neutrino runs. Therefore the inclusion of both resolves the 
degeneracy and hence it is not present in most setups discussed in the literature. For single ion Beta Beams and 
electron capture studies [33-35, 40], this degeneracy needs to be resolved by the experiment or constrained to values 
of 9i3 larger than near future experimental limits. This degeneracy shall hereafter be referred to as the 'energy 
degeneracy' to distinguish it from the usual intrinsic degeneracy which has a different origin. 

Ideally, we would want the energy degeneracy to not be present in the data. To find a condition for this to be so, 
first rewrite the probability (Eq. 3.2) as 

(3.8) 



Pefi — Ii sin^ 26*13 + I2 sin 26*13 cos 5 Ij, sin 26*13 sin 5 + 1^ 



so that all the non-essential constants are tied up in the Ii. As a first step, we attempt to solve Eq. 3.7 for mono- 
energetic neutrinos only with energies Ei and E2. Labelling the respective coefficients as 1} and if we obtain the 
relation 



1 



(sin^ 26li3 - sin^ 26*13) 



1 
I^ 



^3 



(sin 26*13 cos 6*13 cos S — sin 26*13 cos 6*13 cos 5) = 



(3.9) 
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For the energy degeneracy to be resolved, we must have ^13 = ^13 which is true if either 
1. Il/Il = li/ll, i.e. Amli/AEi and Amii/4^2 differ by tt; 
2.ll/li=lf/ll 

In general ^ Ii/^i so we are left with the first condition. The set of equations (3.7) can accommodate the 

solution 

5 = 7r-^-2tan-i (^^^ , (3.10) 

in addition to the trivial solution 6 = 6 for the case ^13 = ^13. Therefore the energy degeneracy for a pair of neutrino 
energies is only completely resolved if condition 1 above and 

^ = ^ — tan~^ 

hold. Therefore, the energy degeneracy is present in general for two beams of mono-energetic neutrinos; however, this 
is not necessarily a nuisance. In section V, it will be found that a combination of mono-energetic neutrino beams 
placed almost on first oscillation maximum and on second oscillation maximum provides some of the largest CP- 
violation sensitivity coverage. In such a case, the CP-even part of the probability vanishes and the energy degeneracy 
is located at 5 ~ n — 5. Since this change will leave the probability invariant, the degenerate region will be the same 
strength and will be symmetrically placed about the J = 7r/2or— ■7r/2 lines. The two regions will always be either be 
CP-conserving or CP-violating at the same time. 

In [35, 40] only neutrinos were used. In those studies, the strategy was to exploit the energy dependence of the 
oscillation signal to break degeneracy and push for a good physics reach; specifically though the combination of bins 
centred around first and second oscillation maximum. The above argument says that the combination of the maxima 
is insufficient to completely break the energy degeneracy, but the 0i3-part of the degeneracy is broken. The reason 
why the degeneracy was only present for very small values of sin^ 2^^i3 in these studies was because the data was 
binned. If one thinks of the data set as predominantly pairs of bins separated by SAm^j^/AE — tt, the location of the 
energy degeneracy is different for each pair and only the true solution is statistically significant. Or more simply, from 
Eq. 3.9, the combination of multiple energies completely breaks the degeneracy as its location is energy dependent. 
In this paper, only two electron capture boosts will be used and so this degeneracy is in general present. 

IV. A DUAL BOOST ELECTRON CAPTURE MACHINE 

In this section, a summary of the simulations carried out will be given. First, the two Lorentz boosts will be 
varied independently to identify pairs of interest. To begin with the number of useful decays per year will be fixed at 
iVions = 10^^, irrespective of the boost, and no atmospheric backgrounds will be included. Since both these parameters 
are critical for the projective CP-sensitivity, they will be introduced shortly once pairs have been identified for further 
study. The neutrino oscillation parameters have been set to the current values [11]. For all simulations, the Gaussian 
form of the function is used: 

j Ci + (/sys ■ Ci)' 

where is the 'true' event rate for boost ji and Ci is the corresponding 'test' event rate. Here /sys = 2 % is an overall 
normalisation error. An intrinsic background of 0.1 % of the unoscillated flux is included to take into account neutral 
current pion production and electrons misidcntified as muons. Inclusion of atmospheric neutrinos will be discussed in 
section VI. Throughout, a 440 kton Water Cerenkov is used with no energy reconstruction. This allows for all charge 
current events to be taken as signal, as oppose to just the quasi-elastic events. The downside of this approach is that 
all beam and atmospheric backgrounds that pass the cuts will be included as signal regardless of their energy. 

The sensitivity to CP-violation of an electron capture machine as a function of two boosts, 71 and 72, and the run 
time fraction, /, will be explored. For a given boost pair (71,72), number of targets jVr, and run time fraction /, the 
simulated event rate is 

POO 

nj, =AfTmt $(-E,) a{u^)P,^^,^ 5{E - E,) dE , (4.2) 
Jo 
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where 771 = /, ?72 = 1 — / and t is the total run time. Although, the CERN-Canfranc baseline is relatively short, the 
full oscillation probability was numerically simulated including matter effects. 

For an electron capture machine, one has complete freedom (within technological bounds) to choose the Lorentz 
boosts of the ion. Clearly, while some choices of boosts will return good sensitivity to CP- violation for given ranges of 
sin^ 2^13, others will result in poor sensitivity. Further, although the boost choice in [34] returned a very competitive 
minimum sin^ 26'i3, the sensitivity region was asymmetric in i5 — >■ —S and there existed degeneracy for sin^ 26*13 ^ 10^^ 
and S < 0. Generating and then preparing CP-sensitivity plots for many pairs (71,72) to compare visually is both 
computationally intensive and a cumbersome undertaking. The task at hand, therefore, is to adopt a measure to select 
boost pairs that return a good overall CP-violation reach or those that are of interest for particular ranges of sin^ 2(?i3. 
One approach would be to determine the minimum sin'^ 2^13 at a given confidence level for a particular experimental 
parameterisation; however, for an electron capture machine, this approach is not desirable for two reasons. As 
mentioned previously, in the previous study [34], the CP-sensitivity region was asymmetric with considerably less 
sensitivity for S > 0. Secondly, there were regions, especially for sin^ 26'i3 ~ 10^^ and 6 < 0°, where degeneracy 
persisted. Choosing a particular value for S then searching for the pair with the minimal sin^ 2^13 will not take into 
account these features. It could be the case that a boost pairing with sensitivity to CP- violation down to very small 
sin^ 2^13 has very poor sensitivity at larger sin^ 2^13 . 

In this article, a measure referred to as the 'CP-coveragcr (i5cov) is adopted. This measure is essentially an ap- 
proximation to an 'integrated CP-fraction', and is described below. For a particular pair (71,72), the CP-coverage is 
determined in the following manner 

1. For sin^2^i3 = S2;i3 € [-s^is, s^is] , CP-violation is tested; for each pair {s2;i3,S) on a grid for the S2;i3 range 
defined previously and S e [—180°, 180°] , a fit to both 6 = 0° and S = 180° is attempted. Solutions for the 
opposite sign(Am3i) are also tested. The minimum (Xmin) calculated is taken. 

2. In this article, 2 degrees of freedom and 99 % confidence levels are adopted. If Xmin > 9-21 for a particular pair 

(^2:13,(5), then 1 is added to a tally. 

3. The fraction of grid points for which CP-conservation can be ruled out is the 'CP-coverage' for chosen range 

S2;13 € [S2;13>S2;13j- 

This procedure is repeated for a grid of (71, 72) pairs. The CP-coverage is a discrete approximation to the integrated 
CP-fraction. Let the CP-fraction, Jjf = <^.f(s2;13) be the fraction of 5 for a given sin^ 2^13 for which CP-conservation 
can be ruled out. Then 

1 /"""is 

(^cov >■ ^ / C^S2;13 (4.3) 

as the number of grid points taken in the range S2;i3 € [s^'il, s^if] and 8 e [—180°, ISO''] tends to infinity. Here, A 
is the area of the region defined by [s2?i3, S2J13] and b e [—180°, 180°] . The CP-coverage is thus an approximation to 
the fractional area of the region in (s2;i3, 5) space for which CP-conservation can be ruled out. 

In the next section, CP-coverage will be used in two stages: firstly for 10"''' < sin^ 26*13 < 10~^; and then for 
each individual decade. Choosing a large range of integration allows one to identify boost pairs that may have the 
desirable CP-features; namely discovery reach down to sin^ 2^13 ~ 10~^ and minimal degeneracy for larger values. 
However, this choice is by no means perfect since if one requires a machine configured for sin^ 2^13 G [lO"'', 10"^], 
say, the information for sin^ 26^13 < 10"^ can distort the relevance of the result. Specifically, it is quite possible that a 
facility with a good minimal sin^ 26^13 may have a poor return with regards CP- violation for large ranges of sin^ 2^13. 
Therefore, to search for further pairs of interest, it is necessary to repeat the analysis for smaller sin^ 2^13 intervals. 

V. RESULTS 

In Fig. 2, the 5 %, 15 %, 25 % and 35 % coverage contours are presented for the case / = 0.3. There are two distinct 
regions in the (71, 72) plane; the first corresponds to the placement of the neutrino flux at the highest energy allowed in 
combination with a flux at second oscillation max;imum. This combination corresponds to the naive expectation that 
good sensitivity should result from combining first and second oscillation maximum. However, there is a (additional) 
large region of (71, 72) space in which sizeable coverage can be achieved. Specifically, for a high boost, 260 < 71 < 400 
and a narrower range of low boosts, 150 < 72 < 185; there is at least 25 % coverage of the (^13, S) plane. This rises 
to 35 % for (71,72) ^ (280,160). This result is repeated for the cases / = 0.4 and / = 0.5 with little variation. 
The small region with 35 % coverage at (280, 160) is the best CP-coverage of the run time fractions and boost pairs 
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FIG. 2: For / = 0.3, the CP-coverage as a function of 71 and 72. 
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FIG. 3: CP-violation sensitivity at 99 % for the pairs (71,72)= (440, 150) (left) and (71,72)= (280, 160) (right). 



studied. The / — 0.3 case runs at 7 = 150 for the longest period, and hence is the easiest technological option. This 
case will be considered henceforth. 

In Fig. 3 the CP-sensitivity plots are shown for the two pairs (440, 150) and (280, 160). These correspond to 
the centre of the 'naive' choice of boosts and maximum coverage region, respectively. It is seen that the first option 
produces an asymmetrical sensitivity region, the best sensitivity for 6 < 0°; any degeneracy is resolved in larger regions 
of parameter space. However, the minimal sin^ 2^13 is larger for S > 0° and this choice of boosts. For S < 0°, the 
lack of sensitivity around sin^ 2^13 ~ lO^'^ that was present in the original study [34] is not present and the minimal 
sin^ 2013 for which CP-conservation can be ruled out at 99 % is also slightly smaller. 

From the right-hand panel of Fig. 3, it is seen that the combination of two lower boosts results in a more symmetrical 
sensitivity region in the (6*13, 6) plane. For large values of sin^ 2^13, the effects of degeneracies are still present, especially 
for S < 0°. This is of little concern, however, since measurement of sin^ 2^13 > 10~^ will be explored with the next 
generation reactor and accelerator long baseline experiments. The minimum sin^ 2^13 for which CP- violation can be 
established is now a factor of 4 larger. This feature is a result of a large region of equivalent solutions being present 
(at 99 % confidence level) at small values of sin^ 26*13 (smaller boosts imply smaller event rates which weaken the 
sensitivity). By choosing two relatively small boosts, the setup has been configured to explore the interference features 
of the appearance probability. In doing so, there is now a poorer resolution on ^13 which comes predominantly from 
the atmospheric features where there is little degeneracy between 6*13 and S. 

In Fig. 4, 90 %, 95 % and 99 % C.L. 2 parameter fits for the pair (440,150) are presented. The overall sensitivity and 
the contributions from each boost are shown and have been computed on the assumption of normal mass ordering. 
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FIG. 4: 90%, 95% and 99% C.L. 2 parameter fits for (71,72) = (440,150) for the CERN-Canfranc baseline (650 km). Plots 
have been produced on the assumption that sin^ 2613 = 10"'^ (left), 10"^ (right), S = 60° (top) and S — —60° (bottom). The 
blues curves correspond to 7 = 440, the red to 7 = 150 and the black to the overall sensitivity. Normal mass ordering has been 
assumed. 
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FIG. 5: 90%, 95% and 99% C.L. 2 parameter fits for (71,72) = (280,160) for the CERN-Canfranc baseline (650 km). Plots 
have been produced on the assumption that sin^ 2^13 = 10~^ (left), 10"'^ (right), 5 = 60° (top) and S = —60° (bottom). The 
blues curves correspond to 7 = 280, the red to 7 = 160 and the black to the overall sensitivity. Normal mass ordering has been 
assumed. 
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FIG. 6: For / = 0.3, the CP-coverage as a function of 71 and 72 for the individual sin^ 2^13 decades. The red circle indicates 
the pair that returns the largest CP-coverage in each analysis. 



The four true value pairs (sin^ 26*13, (5) = (10"^60°), (10"^-60°), (10"^ 60°) and (10^2 -60°) arc shown. The 
choice of boosts corresponds to placing the electron capture flux on second oscillation maximum and, approximately, 
first oscillation maximum. The quantities lS.rn^^L / AEi and lS.rn^^L / AE2 are therefore tt out of phase with each other. 
The sinusoidal shape in solely the cause of interference effects; the solar and atmospheric features of the appearance 
probability have no 5 dependence. In the lower right panel, the energy degeneracy is seen to be explicitly present. 
For the choices 5 = 60° and 6 = —60°, one would expect to see degeneracy at the same ^13 and S = 120° or —120°. 
These degeneracies are in fact present but the large sensitivity regions for the 7 = 150 run have caused them to merge 
into one large sensitivity region. 

The 2-parameter fits for (280,160) are shown in Fig. 5 for the same pairs of true values as before and for normal mass 
ordering. The good CP-sensitivity for low boost pair is a result of the CP-features of the appearance probability being 
out of phase with each other. In particular, the inflexion point of the sensitivity region for one boost corresponds to 
the turning point of the other. There are, in general, two such points; the second corresponds to an energy degeneracy 
at large sin^ 2^13. This solution, however, can be ruled out by near future reactor and accelerator experiments. Note 
that any energy degeneracy that remains does not cross the lines 5= 0° or 5= 180° and therefore does not interfere 
with the sensitivity. The relatively poor sensitivity to CP-violation, is therefore not a consequence of any energy 
degeneracy; instead it can be attributed to poor sin^ 2^13 resolution and low event rates. 

In Fig. 6, the analysis is extended to search for further boost pairs of interest through the narrowing of the 
range of sin^ 2^13 integration. The range sin^ 26113 G [10^"*, 10^^] has been split into three and the scan over the 
(71 J 72) plane repeated for each. It is apparent that the two pairs identified previously are the preferred choices for 
sin2 20i3 € [10-3,10^1]. For sin^ 26*13 € [10""*, 10"^] a new boost pair is identified: (370,170) which is explored 
further in Fig. 7. The CP-violation discovery plot (left panel of Fig. 7) indicates that discovery down to a minimal 
sin^ 2013 = 10~^ can be obtained for both S < 0° and S > 0°. However, this is done at the expense of poor CP- 
violation discovery for all S > 0°. In the right-hand panel of Fig. 7, the CP-discovery analysis is compared to the 
corresponding analysis for (440,150) and (280,160) for the range 10^** < sin^ 26*13 < 10~^, that is to be of interest 
for a Beta Beam-type facility. The reach is good for S < 0° with considerable improvement on the (280,160) pairing. 
It is clear though that the (370,170) is not a desirable option for 5 > 0° as two-thirds of those S, for all sin^ 26*13, 
cannot be distinguished from 5 — 0° or S = 180°. This feature returned a coverage 10 % smaller from the maximum 
(280,160) in the first analysis and hence was not identified there. 
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Hereafter, the analysis will concentrate on the two pairs (280,160) and (440,150). The first is of interest if one 
requires to the electron capture machine to be configured to achieve a broad reach down to sin^ 20i3 ~ few xlO""*, 
whilst the second performs better close to the current limit on sin^ 29i3 and the range to be explored by near-future 
experiments. In the next section, these two facilities will be examined further through the variation of the number of 
useful decays and the level of atmospheric background. 



VI. USEFUL DECAYS AND ATMOSPHERIC BACKGROUNDS 



In the previous section, the scan of the (71, 72) plane was carried out assuming that 10^® useful ion decays per year 
will be available and that the number of atmospheric background events is zero. As discussed in Sec. II, this useful 
decay rate is challenging and would require technology and R&D beyond that required for the standard Beta Beam 
ions. Zero atmospheric backgrounds is not a realistic assumption; with no energy reconstruction nor background 
rejection, every atmospheric event that passes the cuts will be misidentified as an i'^ appearance event. Assuming a 
duty factor of 10~^, i.e. only 0.1 % of the decay ring is filled with ions, there will about 0.03 atmospheric neutrino 
events per kton-year [41]. With a 440 kton detector this amounts to 13.2 events per year. This is a constant that 
will need to be added to the number of events for all count rates in the simulation. The effect of the atmospheric 
background is felt at the sensitivity boundary. In these regions, the number of events is the same order as the 
atmospheric background. The inclusion of the atmospheric background reduces the value of the and pushes the 
boundary inwards. This can be seen trivially with the Gaussian form of a . For true event rate rii, test event rate 

and overall systematic /sys, the Gaussian is given by 

^2^^ {n^-^r ^_ (6.1) 

Inclusion of a constant background, B, is the replacement Ui n.i + B and similarly for Q. This constant translation 
of the event rates leaves the numerator of the unchanged whilst increasing the denomenator. The result is that 
the is reduced. 

The physics reach of the electron capture beam cannot be stated until the atmospheric background is included. 
Further, it is not known what is a realistic useful decay rate per year - it is reasonable to assume that 10^^ may well 
be close to an upper limit though. A technique to improve the decay rate is to increase the duty factor which in turn 
introduces a larger atmospheric background rate. Since the useful decay rate is not known, the approach requested 
by the experimentalists [42] is to treat the useful decay rate and number of atmospheric events as independent. 

This analysis has been carried out using the GP-coverage measure and sin^ 26'i3 £ [10~^, 10""*^], as previous. Note, 
however, that its use here is for a different purpose. Whereas before the goal was to identify boost pairs of interest, here 
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FIG. 9: CP-violation sensitivity at 99 % confidence level for / = 0.3 and A'atm ~ 13.2 per year for (71, 72) = (440,150) (left) and 
(280,160) (right). In both cases, the sensitivity contour has been constructed for a range of useful decay rates. 



it is used to monitor the area change of the 99 % confidence region. As the decay rate drops and/or the atmospheric 
neutrino level is increased, this region will shrink by shifting inwards. The CP-coverage allows one to adequately 
measure this change. 

The results are presented in Fig. 8. The two pairs (440,150) and (280,160) with / = 0.3 have been re-simulated 
incorporating the number of background events per year and varying the number of useful ion decays. It is seen 
that the electron capture machine, like all Beta Beam type machines, is statistics dominated; the coverage of the 
(sin^20i3, 5) plane is reduced considerably with drops in the decay rate. For example, if the average number of 
atmospheric background events per year is 10, reducing the useful decay rate by a factor of 2 reduces the coverage by 
almost 10 % of the plane for both cases. When the coverage is only 35 % for 10^^ useful decays, this is a substantial 
drop. On the other hand, if one instead opts to regain that factor of 2 by loosening the duty factor, i.e. double the 
number of atmospheric neutrino events, the reduction is only 2-3 %. This observation seems to agree with a previous 
study for a CERN-Frejus ^^Ne and ^Be Beta Beam [43] where it was demonstrated that a loosening of the duty factor 
does not lead to the substantial reductions in physics return as one might assume. 

To see the results more explicitly, in Fig. 9 CP-sensitivity plots have been constructed (for both cases) including a 
constant atmospheric background of 13.2 events per year for a range of useful ion decay rates. The choice of boost 
pairs presented here provide excellent resolution of the energy degeneracy. As a consequence, a lower event rate merely 
reduces to overall sensitivity of the setup, rather than let degeneracies ruin the sensitivity. The first manifestation of 
any degeneracy is for TVions = 2 x lO^'^ per year for (440,150) at sin^ 2^13 ^ 10^^ and S< 0° . In Fig. 10 this exercise 
has been repeated but with Anions = 10^^ per year kept fixed and the atmospheric event rate varied between and 
100 events per year. Here the reduction in sensitivity is less extreme with no degeneracy appearing. The effect of 
the background is only felt through the denomenator of the so the facilities ability to resolve degeneracy is largely 
unaffected. An increased background therefore pushes the sensitivity boundary at a given confidence level inwards 
towards the set of pairs (6'i3, 5) where the event rate and background rate are of the same order. 
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FIG. 10: CP-violation sensitivity at 99 % confidence level for / = 0.3 and Anions = 10^* per year for (71, 72) = (440,150) (left) 
and (280,160) (right). In both cases, the sensitivity contour has been constructed for a range of atmospheric neutrino events. 

These plots indicate that the physics reach of the electron capture machine is a struggle between sourcing a plentiful 
useful decay rate whilst keeping the atmospheric background to a minimum. With binning of the neutrino signal, 
the constant atmospheric neutrino rate can be reduced whilst keeping the useful decay rate fixed. There are two 
drawbacks to this possibility: 

1. The use of QE-events would be necessary 

2. An event reconstruction efficiency needs to be included 

The electron capture machines need to use large Water Cerenkov detectors to achieve a competitive event rate: the 
cross-sections are small at the lower energies studied at short baselines. The energy reconstruction of neutrino events 
at a Water Cerenkov is through selection of quasi-elastic events with an efficiency of about 60-70 %. This is in addition 
to the loss of events from not using all the charge-current events. There will also be a small further loss owing to finite 
energy resolution. Binning the signal will reduce background but the reduction in events is essentially equivalent to 
using a smaller useful decay rate. There is therefore likely to be no gain in using this approach. 

Alternative strategies to reduce the atmospheric background rate would be to either introduce cuts on the signal or 
to use directional information. Since the neutrinos are at a known energy, one can impose cuts on the signal for larger 
and smaller energies. The atmospheric events are skewed towards low energies, so this would work best for large ion 
boosts. At high energies, the neutrino direction is correlated with the observed lepton direction. This gives a means 
to separate the beam signal from the approximately isotropic atmospheric signal. However, any gain from these extra 
reductions will probably be offset by an increase in the duty factor in an attempt to achieve a useful decay rate. 

VII. MONO-ENERGETIC ANTI-NEUTRINOS 

The electron capture machine setups with the largest CP-violation coverage are either asymmetrical in 5 or sym- 
metrical, but with a larger sin^ 2^13 limit. It is clear that unless vast improvements in the available number of useful 
decays per year are achieved or envisaged, the physics reach of electron capture machines cannot compete with the 
standard Beta Beams over equivalent baselines which can rule out CP-conservation down to sin^ 26*13 ^ lO"'' and are 
close to symmetric in d. The reason for this is the use of a neutrino run only. The need for information complemen- 
tary to first oscillation maximum data, forcing runs at smaller boosts, has reduced the overall event rate across the 
whole experimental run. One option to explore is whether it is possible to construct a facility with a mono-energetic 
anti-neutrino flux at flrst oscillation in combination with an equivalent neutrino flux from electron capture. The two 
fluxes would be complementary with the need for additional runs at lower energies unnecessary. The event rate across 
the entire run of any experiment would be higher, with less susceptibility to variations in the useful decay rate. One 
process put forward as a complementary source of mono-energetic anti-neutrinos is the bound beta decay process 
(BBD) [47]. 

In a continuum beta decay (CBD) the electron is virtually always ejected from the atom with capture in the outer 
orbitals strongly suppressed due to weak bindings and small wavefunction overlaps. However, for a fully ionised atom, 
electron capture in the K-shell (the bound beta decay process (BBD)) can have a substantial branching ratio, provided 
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the transition Q-value is sufficiently small and proton number of the ion is large: 




Ql\MO)\' 
f{Qc,z) 



(7.1) 



where 




meiQc- E fF{Z, E)dE 



(7.2) 



is the integral over phase space for CBD and V'n(O) is the wavefunction of the electron in the nth orbital. For a fully 
ionised atom, the Q-value for the BBD process is given by 



where |Ai3*°*(/', /)| is the difference in binding energies of the complete parent and daughter atoms, and |i?„(/')| is 
the binding energy of the electron captured into the orbital n of the daughter nucleus. Specifically, the Q- values for 
the BBD and CBD channels are separated by |B„(7)| in the rest frame of the nucleus. 

In the last couple of years, use of BBD has been forwarded for very short baseline experiments [44] searching for 
oscillations of so-called 'Mossbauer neutrinos' [45]. The possible reaction rates for a range of candidate ions sourced 
from a low-boost Beta Beam has also be discussed in [46] . Here the possibility of using BBD as a neutrino source for 
a long baseline neutrino oscillation experiment is considered as was suggested in [47] where the authors proposed the 
use of ions that can BBD as well as electron capture and CBD. The paper defined the 'CP-evenness' as 



where T are the fluxes and a the cross-sections at definite energies. Since T{i>)a{i>) is the un-oscillated number of 
events at the detector, it was suggested that an optimum neutrino beam is one with 77 « (equal un-oscillated neutrino 
ans anti-neutrino events). The feasibility of this approach and the combination of an electron capture machine and 
bound beta decay machine was discussed in [48]. Below I will summarise the main reasons why a machine based on 
the BBD process is not a realistic proposition. 

1. The branching ratios are too small. For a fully ionised ion to have a BBD substantial branching ratio requires a 

large proton number and a very small Q-value. There is always a substantial contribution from the CBD. A scan 
of a database searching for suitable ions with half- lives in the range 0.5 sec < < 8 min and a dominant decay 
channel returned very little. The largest BBD branching ratio found is for 207rpj8H-. ^2 %. Therefore to source 
an equivalent anti-neutrino event rate requires an extra factor of 10 useful decays, compared to electron capture 
decays, and a further factor of 3 to compensate for the lower anti-neutrino cross-section. The ions suggested for 
electron capture and BBD in [47] had much lower branching ratios around 1 %. 

2. Magnetised detectors are necessary for ions with both electron capture and BBD. For ions that electron capture 
and BBD, it is necessary to separate fi~ and /i"*" events at the detector. For MIND detectors [49], the threshold 
is too high for the intermediate baselines of Europe. Magnetised Liquid Argon or Totally Active Scintillator 
detectors would be necessary. However, the previous comment indicates that a factor of 100 is missing to get a 
useful number of mono-energetic neutrinos of order 10^*, so this is a mute point. 

3. Very large boosts are required to separate BBD and CBD channels in the laboratory frame. CBD will always 
be the dominant channel. It is therefore necessary to separate this out from the BBD channel otherwise the 
mono-energetic nature of the BBD is not being exploited. Failure to do this would result in a Beta Beam using 
a high proton number ion with only the endpoint of the flux in use. The energy split between CBD and BBD 
in the laboratory frame is E = 2jBi{Y), for boost 7 and decay to the K-shell of daughter ion Y. For example, 
207rpj8i-i- ^g^g Bi{Y) = 99 keV. For a boost of 7 = 400, this corresponds to a required binning and energy 
resolution AE < 80 MeV, which is tough. More conservative binnings of 150-200 MeV and greater are often 
considered. Boosts in excess of 800 will be required to have an empty bin between the CBD and BBD fluxes 
in these cases. Such fluxes require the LHC and will diminish the useful decay rate further due to its leading 
order I/7 dependence. For all ions with smaller proton numbers, the experimental task will be harder owing to 
the smaller binding energies. 



Q%+ = Qc +1B„(/')|-|AB*°*(/',7)1 
= Qg+ + lB„(/')| , 



(7.3) 



T{ue;E)a{u^-E) - T{v^; E)a{V^- E) 
T{ve;E)a{Vf,-E)+F{ve;E)a{vf,;E) ' 



(7.4) 



15 



In short, the implausibihty of achieving 0(10"'^^) useful mono-energetic anti-neutrinos per year and separating out the 
BBD anti-neutrinos from the CBD anti-neutrinos does not make a Beta Beam-type facility based on BBD a viable 
future facility. This does not mean that anti-neutrinos cannot accompany an electron capture beam. That can be 
achieved by combining an electron capture machine with the standard anti-neutrino ions; ^He and ^Li. This is left 
for a different study. 

VIII. SUMMARY AND CONCLUSIONS 

A electron capture beam provides a means to source long baseline neutrino oscillation experiments with a mono- 
energetic Ve flux. By altering the Lorentz boost of the ion, one can freely choose the laboratory neutrino energy up to 
some maximum determined by the Q- value and the maximum energy of the accelerator. Since a single neutrino energy 
is insufficient to determine the unknown neutrino mixing parameters and resolve any degeneracy, it is necessary to 
include at least two boosts in any physics strategy. In the present article, a previous work [34] examining the physics 
reach of an electron capture beam sourccd from CERN and directed towards a 440 kton fiducial mass Water Cerenkov 
detector was extended to include the study of the choice of boost, the relative run times of each boost, the number 
of useful ion decays, and the number of atmospheric neutrino events. The use of neutrinos only, but with multiple 
energies introduces a variant on the intrinsic degeneracy that is not present in studies that use both neutrino and 
anti- neutrinos. This 'energy degeneracy' was discussed and shown to present in electron capture beam facilities unless 
the two oscillation phases Am|]^/4i<^ arc separated by tt and the CP-phase 5 takes a specific value. 

In the first instance, the number of useful decays was fixed to 10^^ per year, and atmospheric backgrounds were 
neglected. Little variation was found for the relative run times; however, two regions of the (71 , 72 ) plane returned large 
CP-coverage. The naive pairing of (almost) first and second oscillation maximum had a CP-coverage of approximately 
35 % but was asymmetric in 5 — >• —8. For (5 < 0, this choice could rule out CP-conservation down to sin^ 2^13 = 10~^, 
but only sin^ 2^i3~ 2 • 10^^ for 5 > 0. There exists a slightly larger coverage for the pair (280, 160) with run time 
fraction / = 0.3. This choice returns a (roughly) symmetrical CP-sensitivity region with CP-conservation ruled out 
down to sin^ 26*13- 4 • 10"^. 

No R&D has been carried out on the feasibility of sourcing lO""^^ useful decays per year for high proton number ions. 
A brief study adapting Beta Beam codes indicated that the decays are too slow and large tune-shifts would need to 
be accommodated. With existing technology, the number of useful decays is two orders of magnitude too small. At 
best, 10^^ useful decays per year appears to be overly optimistic and should be considered as the hard limit on the 
yearly rate. Some of this problem could be alleviated by keeping large numbers of electrons bound to the nucleus 
whilst still having a large enough charge-to-mass ratio to reach the sought boost. 

In a standard Beta Beam, a restrictive duty factor of Sf — 10~^ is normally taken as necessary to reduce the 
atmospheric neutrino flux to a level that does not destroy the sensitivity. The next step, therefore, was to examine 
the CP-sensitivity as a function of both the yearly useful decay rate and the number of atmospheric neutrino events 
in the appearance sample. In both cases, the results were more volatile to small changes in the useful decay rate than 
the atmospheric background. For instance, taking iVjons = 10^^ and i?atm = 10 events per year, the CP-coverage is 
approximately 35 % in both cases. If the decay rate is dropped by a factor of 2. then the coverage falls to about 
23 %. If instead the atmospheric background is doubled, then the coverage only drops to about 32 %. The principle 
reason for this difference is the total event has manifestations in both the numerator and denomenator of the 
function whereas the background only affects the denomenator. The discrepancy between two (Q\-j,^S) is felt though 
the difference squared in the numerator. Drops in the total event rate can therefore lead to degenerate solutions 
dropping below the required statistical significance of the test thus reducing the sensitivity in a manner beyond 
simply a scaling of the . The background merely reduces the without interfering with its ability to rule out 
degenerate solutions. Its manifestation is simply to push the sensitivity contours inwards. 

Ultimately, the limiting factor for an electron capture machine is the number of useful decays. The advantages of 
using a mono-energetic neutrino beam is the precise knowledge of the neutrino energy and ability to choose which 
part of the oscillatory structure of the appearance probability is to be explored. Ironically these are also the main 
disadvantages. Since a single energy is insufficient to extract the unknown mixing parameters, it is necessary to 
run for substantial periods of time at low energies where the event rate is small. The availability of an equivalent 
anti-neutrino beam would remove this necessity. Ions that bound beta decay could source such a beam; however, this 
will not be practical with the available technology. 

The feasibility of an electron capture beam is an open question. However, the most optimistic parameterisation 
considered in this article is likely to be beyond the hard limits imposed by available technology. It must be conceded 
that if 013 were to be very small, the electron capture will not be a competitive facility, unless a fast decaying ion that 
can be produced in large quantities is found. This can be seen explicitly in Fig. 11 where the CP-discovery for the 
boost pairs (280,160) and (440,150), and 20 atmospheric events per year, are compared to the 7 = 350,350 Beta Beam 
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FIG. 11: CP-violation sensitivity at 99 % confidence level for / = 0.3 and A'atm = 20 per year for (71, 72)=(440,150) (black 
solid) and (280,160) (blue dot-dashed). The red dashed line displays the physics return for a standard Beta Beam directed 
along the CERN-Canfranc baseline with a boost 7 — 350 for both neutrinos and anti-neutrinos. 

introduced in [41] with the baseline shortened to CERN-Canfranc, as in this article. The simulation of this facility 
uses 1.1 X 10^* neutrinos per year sourced from the decay of ^^Ne, and 2.9 x 10^^ anti- neutrinos per year sourced 
from decays of ^He. The run time is 10 years with an equal split between neutrino and anti-neutrino running. We 
see that for a (440,150) electron capture facility, the physics return is poorer but relatively competitive for S < 0°. 
For these 6, the (280,160) facility is only competitive for sin^ 26*13 e [10"3,10-2]. For S > 0°, the (440,150) facility 
is outperformed by the standard Beta Beam by 2 orders of magnitude whilst, again the (280,160) facility is only 
competitive for sin^ 26*13 € [10"^, 10~^]. 

In conclusion, an electron capture machine is unlikely to be a competitive facility for the measurement of CP- 
violation. Taking into account technological challenges and restrictions, the CP-violation reach of a possible electron 
capture machine is probably going to be limited to sin^ 26*13 > 10~^. This is not to say that an electron capture 
machine is a redundant option. The sensitivity suffered from the need to run for substantial periods of time at low 
energies. If the electron capture machine is used in combination with another facility, this will be become unnecessary. 
In such a scenario the role of the machine would be to complement a facility that has both neutrino and anti-neutrino 
fluxes and/or coverage of a large energy range. In this case, the electron capture machine could be run for 100 % of 
the time at a well chosen energy to compensate for any deficiency in the physics reach of the other facility. 
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